Sudden Hadronization in Relativistic Nuclear Collisions 
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We formulate and study a mechanical instability criterion for sudden hadronization of dense matter 
fireballs formed in 158^4 GeV Pb-Pb collisions. Considering properties of quark-gluon matter and 
hadron gas we obtain the phase boundary between these two phases and demonstrate that the 
required deep QGP supercooling prior to sudden hadronization has occurred. 
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Hot and dense hadron matter fireball is formed in cen- 
tral collisions of relativistic 158^4GeV heavy Pb- nuclei 
with Pb-target, comprising a new state of matter [Q. 
Driven by internal pressure, a fireball expands and ul- 
timately a breakup (hadronization) into final state par- 
ticles occurs. Early on in the hadron production data 
analysis it was discovered that strange hadrons emerge 
from a source in which strange s and antistrange s quarks 
have same size phase space [0. This is the case if final 
state hadrons are directly produced by the deconfmcd 
quark matter phase. 

The required sudden fireball breakup could arise if 
a fireball made of a the new form of matter signifi- 
cantly supercools, and in this state encounters a strong 
mechanical instability ||. Despite extensive ensuing 
study, the mechanisms determining when and how sud- 
den hadronic particle abundances are formed (chemical 
freeze-out) have not been fully understood 0-k|. How- 
ever, there is growing evidence that the fireball breakup 
occurs over a relatively short period of time |t] ^] . 

We propose and study here a natural mechanical in- 
stability criterion ensuing the fireball expansion into the 
metastable supercooled state. We consider the exploding 
fireball dynamics in its center of momentum frame of ref- 
erence. The surface normal vector of exploding fireball 
is n, and the local velocity of matter flow v c . The rate of 
momentum flow vector V at the surface is obtained from 
the energy-stress tensor T^i p0| : 



P = P«n+(P«+£«)^^ 

1 — v/ 



(1) 



The upper index (?) refers for the intrinsic energy density 
e and pressure P of matter in the frame of reference, 
locally at rest, i.e. observed by a co-moving observer. 
We omit the superscript (i) in the following. For the 
fireball expansion to continue, V = \P\ > is required. 
For V — > at v c ^ 0, we have a conflict between the 
desire of the motion to stop or even reverse, and the 
continued inertial expansion. 

When the flow velocity remains large but V — > 0, the 
intrinsic pressure P must be negative. As illustration 
consider the fireball to be made of a quark-gluon liquid 



confined by an external vacuum pressure B. The total 
pressure and energy comprise particle (subscript p) and 
the vacuum properties: 



P = P P -B, e = e p + B. 
Eq. (0) with V = thus reads: 



Bn = P p n 



v c v c -n 
' 1 - vl 



(2) 



(3) 



and it describes the (equilibrium) condition where the 
pressure of the expanding quark-gluon fluid is just bal- 
anced by the external vacuum pressure. 

Expansion beyond V — > is in general not possible. 
A surface region of the fireball that reached it but con- 
tinues to flow outwards must be torn apart. This is a 
collective instability and thus the ensuing disintegration 
of the fireball matter will be very rapid, provided that 
much of the surface reaches this condition. We adopt the 
condition V = at any surface region to be the instability 
condition of an expanding hadron matter fireball. 

Negative internal pressure P < is a requirement. At 
this stage the fireball must thus be significantly super- 
cooled. The adiabatic transfer of internal heat into ac- 
celerating flow of matter provides the mechanism which 
leads on the scale of r = 2 10~ 23 s to the development of 
this 'deep' supercooling 

It is possible to determine experimentally if the condi- 
tion P < has been reached. Namely, the Gibbs-Duham 
relation for a unit volume: 



P = To + n h is h - e , 



(4) 



relates the pressure, to entropy density a = S/V , energy 
density e = E/V , and baryon density = b/V , V is the 
volume, T is the temperature, and /ib the baryochemical 
potential. Dividing by e we obtain: 



PV = J\_ t 
E E/S E/b 



- 1 . 



(5) 



The microscopic processes governing the fireball breakup 
determine how the quantities entering the right hand side 
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of Eq. (||) are changed as hadrons emerge. Understanding 
this we can determine, if the intrinsic fireball pressure 
prior to breakup, has been negative. 

The energy E and baryon content b of the fireball are 
conserved. Entropy S is conserved when the gluon con- 
tent of a QGP fireball is transformed into quark pairs 
in the entropy conserving process G + G — > q + q. 
Similarly, when quarks and antiquarks recombine into 
hadrons, entropy is conserved in the range of parame- 
ters of interest here. Thus also E/b and S/b is conserved 
across hadronization condition. The sudden hadroniza- 
tion process also maintains the temperature T and baryo- 
chemical potential /i;, across the phase boundary. What 
changes are the chemical occupancy parameters. As glu- 
ons convert into quark pairs and hadrons "f g — > but the 
number occupancy of light valance quark pairs increases 
Iq > Iqa — 1 increases significantly, along with the num- 
ber occupancy of strange quark pairs 7 S > 7 SD ~ 1. 

The sudden hadronization picture differs from, e.g. , 
the droplet-driven reequilibration transformation C3,ff3l , 
in which chemical equilibrium of valance quark pair abun- 
dances is maintained. Instead a change (reheating) of sta- 
tistical parameters T, ^ occurs, along with a possible for- 
mation of a mixed phase required for volume expansion. 
In such reequilibration hadronization picture, entropy in- 
crease can also occur To draw the line between 
these two hadronization pictures (equilibrium/sudden), 
we need to determine the quark pair occupancy param- 
eters 7j. Our study of final state hadron abundances 
strongly favors ja ~ exp(m^/2T) > 1,7, ~ 7, |fLl|. 
Evaluating Eq. (ph using the results of our data analy- 
sis, we indeed obtain Pf < 0. The magnitude of |P/| 
can vary between a few percent (in terms of energy den- 
sity Ef /V), up to 20% for the latest published result jjl) . 
The precise value, which arises from several cancellations 
of larger numbers is sensitive to the strategy of how the 
currently available experimental data is described, e.g., 
if strangeness conservation is implemented, and if so, if 
differentially at each rapidity, or as an overall conserva- 
tion law; how many high mass resonances can be excited 
in hadronization process, etc ... 

Importantly, we have not been able to obtain a scheme 
of hadron production analysis which describes the data 
with '% 2 /dof < 1 and would not imply P < for the 
hadronizing fireball matter. On the other hand, if we do 
force the hadronizing particles to be in chemical equilib- 
rium, we find x 2 /dof > 2.5, dof = 10 in our analysis 
which agrees for this limit with fl5| , and in this case we 
find P > 0. 

Understanding in detail the breakup condition V — * 
requires that we model the shape and direction of flow 
in the late stage of fireball evolution, obviously not an 
easy task. However, considering n-V — ► 0, we find the 
constraint: 



-PV/E 



(0) 



vectors v c and ft are everywhere parallel, thus k — > 1. 
However, in 158A GeV Pb-Pb reactions the longitudinal 
flow is considerably greater than the transverse flow jl6) , 
and we note k — > for a longitudinally evolving cylin- 
drical fireball. For the Pb-Pb collisions considered here, 
our analysis suggest 0.1 < n < 0.6. 

We now substitute, in Eq. (||), the fireball matter prop- 
erties employing the Gibbs-Duham relation, Eq. (Q) , and 
arrive at: 



§=|T„ 



Mb 

S/b 



1 



(7) 



Eq. (0) establishes a general constraint characterizing the 
fireball breakup condition. 

The solid line, in figure [l], shows the behavior of v c (Th) 
constraint arising from Eq. (^) for the example E/S — 
0.184± 0.05 GeV (error range shown by dotted lines), k — 
0.6. Outside of the region bounded by the solid line (i.e., 
for greater Th and v c ), the flow expansion can occur as 
the internal particle pressure is greater than the confining 
pressure. Also shown in figure [j] is hadron production 
analysis result jllj and its statistical error, the systematic 
error is of same magnitude. The agreement of theory and 
experiment results from the choice of non-spherical flow 
with specific freeze-out shape described by the average 
value k = 0.6, see Eq. (^|). Figure |] illustrates the great 
sensitivity to the analysis on the freeze-out constraint. 
The dashed horizontal line, in figure [l], is the velocity of 
sound of the interacting quark-gluon liquid, which barely 
differs from 1/^/3 @. 
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FIG. 1. Fireball velocity as function of breakup temper- 
ature constraint for the case E/S = 0.185 ± 0.005 GeV, with 
S/b — 42 and fj,^ — 0.2 GeV; the dotted lines describe the 
uncertainty in the determination of E/S. Dashed line: veloc- 
ity of sound of relativistic quark-gluon liquid. Also shown is 
hadron production analysis result fu| . 

We have so far not used in the discussion any key spe- 
cific property of the equations of state of the matter filling 
the fireball. However, our results imply that the matter 
inside the fireball is deeply supercooled. Can this be 
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the deeply supercooled liquid of quarks and gluons? In 
fact a study of QGP equations of state employing prop- 
erties of QCD interactions and thermal QCD JL7J, fine 
tuned to agree with the properties of lattice QCD results 
jl8| suggest that. We extend this study to consider the 
phase boundary. The thin solid line in the T, /ib plane 
in figure ^ shows where the pressure of the quark-gluon 
liquid equals the equilibrated hadron gas pressure. The 
hadron gas behavior is obtained evaluating and summing 
the contributions of all known hadronic resonances con- 
sidered to be point particles. When we allow for finite 
volume of hadrons |l9|] , we find that the hadron pressure 
is slightly reduced, leading to some (5 MeV) reduction in 
the equilibrium transition temperature, as is shown by 
the dashed line in figure ||. For vanishing baryo-chemical 
potential, we note in figure ^| that the equilibrium phase 
transition temperature is T pt ~ 172 MeV, and when fi- 
nite hadron size is allowed, Tf p ~ 166 MeV, The scale in 
temperature we discuss is result of comparison with lat- 
tice gauge results. Within the lattice calculations [pl|, it 
arises from the comparison with the string tension. 
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FIG. 2. Thin solid and dashed lines: equilibrium phase 
transition from hadron gas to QGP liquid without and with 
excluded volume correction, respectively. Dotted: breakup 
condition at shape parameter k = 0.6, for expansion veloc- 
ity vl = 0,1/10,1/6,1/5,1/4 and 1/3, and thick line for 
v c = 0.54. The experimental point denotes chemical nonequi- 
librium freeze-out analysis result pd| . 

The dotted lines, in figure |^, correspond to the condi- 
tion Eq, (|3|) using the shape parameter k = 0.6, Eq. (0), 
for (from right to left) v 2 c = 0, 1/10, 1/6, 1/5, 1/4 and 1/3. 
The last dotted line corresponds thus to an expansion 
flow with the velocity of sound of relativistic noninteract- 
ing massless gas. The thick solid line corresponds to an 
expansion with v c — 0.54. The hadron analysis result is 
also shown O]. Comparing in figure ||thin solid/dashed 



with the thick line, we recognize the deep supercooling 
as required for the explosive fireball disintegration. The 
super-cooled zero pressure P = QGP temperature is 
at T sc = 157 MeV, (see the intercept of the first dashed 
line to the right in figure |^) and an expanding fireball can 
deeply super-cool to T& sc ~ 147 MeV (see the intercept of 
thick solid line) before the mechanical instability occurs. 

Deep supercooling requires a first order phase transi- 
tion, and this in turn implies presence of latent heat B. 
Physical consistency then requires presence of external 
(negative) vacuum pressure —B. More precisely, the vac- 
uum contribution to the physical properties of deconfmcd 
matter can be derived from In Z vac = —BY (5 : 



-B, 



_ T 

*vac — ~y In -^vac 

_ d In Z vac _ 
£vac - vdp - a 



dlnB/B a 

91n/3/A> 



(8) 
(9) 



The temperature, T = 1/(3, dependence of the vacuum 
pressure has been considered within the model of color- 
magnetic vacuum structure |2^,|2l[. Near to the phase 
transformation condition, the variation of B with (3 is 
minimal (see figure 2 in [pfj|| ), and thus the logarithmi- 
cally small last term in Eq. (^|) can be, in principle, ig- 
nored. 

We now combine the theoretical properties of the QGP 
equations of state with the dynamical fireball properties 
in order to constrain B. Reviewing Eq. (jHI), we obtain: 



PV 
~E 



-£qgp 



Pv = B. 



(10) 



To evaluate B, we note that lattice results for £qgp are 
well represented by £qgp = aT A , with a ~ 11, value 
extrapolated for the number of light quark flavors being 
rif = 2.5 at the hadronization point p7[ . We obtain, for 
the fireball formed in Pb-Pb reactions, 

0.2 • HTjf ~ 0.17 GeV/fm 3 < B. 

Is our picture of fireball evolution compelling? We 
found that particle production occurred at condition 
of negative pressure expected in a deeply supercooled 
state and have shown internal consistency with (strange) 
hadron production analysis involving chemical non- 
equilibrium. Moreover, these chemical freeze-out con- 
ditions agree with thermal analysis Q , allowing the con- 
jecture that the explosive quark-gluon fireball breakup 
forms final state hadrons, which do not undergo fur- 
ther reequilibration. However, we noted that the chem- 
ical equilibrium reaction picture differs from ours only 
in terms of its statistical significance (x 2 /dof > 2.5, 
p"5[). It produces a chemical freeze-out temperature of 
T = 168 MeV, just the value we found for an equilibrium 
phase transition implicit in the assumption of chemical 
equilibrium. The higher chemical freeze-out temperature 
produces greater population of excited hadronic states. 
Their decays deform hadron spectra, and this allows for 
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a second evolution stage with a thermal freeze-out tem- 
perature at or below 120 MeV. 

Strongly in favor of the here described sudden QGP 
hadronization resulting in chemical nonequilibrium reac- 
tion picture is the presence of a hadron multiplicity ex- 
cess, related to an entropy excess [p2[ . This is seen both 
as multiplicity per baryon, and as increase of multiplicity 
comparing pA to AA reactions. We could not describe 
this effect by admitting other physical models such as, in 
medium, change of hadron masses. We have found that 
invariably the statistical significance of the analysis de- 
creases as we modify individual hadron properties in an 
ad-hoc fashion, while maintaining chemical equilibrium 
of hadron abundances. Our experience shows that the 
only theoretical description of hadron production data 
that works (% 2 /dof < 1) requires excess of valence quark 
pair abundance, irrespective of the detailed strategy of 
data analysis. This agrees well with the dynamical study 
of nuclear collisions within the UrQMD model concludes 
that at CERN energies the chemical non-equilibrium is 
required to characterize the numerical results in terms of 
a statistical model pj| . 

Also in favor of our result is the conclusion of Csorgo 
and Csernai || : who required as verification for the pres- 
ence of a deeply supercooled state of matter and sud- 
den hadronization: i) short duration and relatively short 
mean proper-time of particle emission, now seen in par- 
ticle correlations J7j, ii) clean strangeness signal of QGP 
p4j ; iii) universality of produced particle spectra which 
are the remarkable features of strange particle production 
p5[ v) no mass shift of the phi-meson; despite extensive 
search such a shift has not been found by the NA49 col- 
laboration | p6[ . 

In summary, we have introduced a constraint, Eq. (0), 
which relates the physical and statistical properties of 
the hadronic fireball at the point of sudden breakup. We 
obtained this result from mechanical stability considera- 
tion, employing only properties of the energy-stress ten- 
sor of matter, and the Gibbs-Duham relation, Eq. 
We showed that this constraint is consistent with analy- 
sis results obtained considering the experimental particle 
production data for Pb-Pb collisions at 158 A GeV. We 
further studied the behavior of the phase transition be- 
tween hadron gas and quark-gluon liquid, and have deter- 
mined the magnitude of the deep supercooling occurring 
in the fireball expansion. Employing a lattice-QCD based 
estimate on number of degrees of freedom in the energy 
density of the QCD thermal matter, we obtained a con- 
straint on the magnitude of latent heat /vacuum pressure 
B > 0.17GeV/fm . We conclude that both in theoreti- 
cal study of the data as well as for reasons of principle 
the deciding factor about the sudden nature of the phase 
transformation is the absence of chemical equilibrium. 
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